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ABSTRACT 


The  MGR[v]  algorithm  of  Ries,  Trottenberg  and  Winter  with  v  =  0 
and  the  Algorithm  2.1  of  Braess  are  essentially  the  same  multigrid 
algorithm  for  the  discrete  Poisson  equation:  -AhU  =  f  .  In  this  report 

CJL 

we  consider^the  extension  to  ther  general  diffusion  equation.  -V  •  pVu  =  f  , 
P  =  p(x,y)  >  Pn  >  0  In  particular,  for  the  two-grid  scheme  we  reobtain 

^  , - - -  .  it  t~  C  ,  h  t  k  •-  •  *  • 


the  basic  result '  p  <  j  (1+Kh)  in  the  stronger  form  ^  <-^0+Kh).r 

Computational  results  indicate  that  other  constant  coefficient  results 


carry  over  as  well.  $  JAAs^r- 
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I.  Introduction 


- )  Multi grid  methods  are  proving  themselves  as  (very)  successful  tools 

for  the  solution  of  the  algebraic  equations  associated  with  discretization 

of  Elliptic  Boundary- Value  problems  -  see  [1],  [3],  [4],  [5],  [9]\  Never- 

v  y**' 

theless,  it  seems  we  are  just  beginning  to  understand^  this  powerful  idea. 
Hence,  there  is  a  need  for  continued  probing,  experimentation  and  new 
proofs  -  less  for  the  sake  of  proof  and  more  for  the  sake  of  insight.  — ) 
In  [2]  Braess  proposed  and  analyzed  a  class  of  multigrid  methods.  In 
particular,  he  considered  a  particular  algorithm  for  the  Poisson  Equation  - 
"Algorithm  2.1".  He  shows  that  the  contraction  number  p  for  a  two-grid 
method  is  given  by 

(1)  P<^! 

This  result  holds  whenever  ft  is  a  polygonal  domain  whose  sides  have  slope 
+1,0  or  «  and  the  discretization  satisfies  an  additional  condition 
(see  ftl  of  section  2).  In  [8]  Ries,  Trottenberg  and  Winter  discuss  the 
class  of  MGR[v]  methods  for  the  Poisson  Equation  in  a  square.  Using  Fourier 
Analysis  they  obtain  an  explicit  formula  for  the  corresponding  contraction 
numbers  p[v].  In  particular,  they  obtain  -  for  two  grids 

(2)  PM-J.  pM-£.  »M-Z^$5rr- 

As  it  happens  MGR[0]  is  the  same  as  the  "Algorithm  2.1"  and  the  results  of 
[2]  and  [8]  are  consistent.  The  results  of  [8]  are  more  precise  for  more 
restricted  problems. 
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In  this  report  we  consider  the  more  general  diffusion  equation 

(3)  -V*p(x,y)Vu  3  f  in  ft, 

u  =  0  on  3ft 
p(x,y)  >  pQ  >  0  and 

where  ft  may  be  a  general  bounded  piecewise  smooth  domain,  or,  ft  is  a 
polygonal  domain  whose  sides  have  slope  +  1,  0,  or  «.  We  employ  the  usual 
five-point  difference  analog  of  (3)  and  seek  to  solve  the  (large)  system  of 
linear  algebraic  equations.  We  consider  a  class  of  linear  multi-grid  methods 
which  include  the  MGR[v]  methods  when  p(x,y)  =  1  .  Our  basic  result  is  the 
following:  Consider  the  two-grid  method.  Then 

P<  ||  MG  ||,  <  i(l+Kh) 

Lh  L 

where  ||  ||.  denotes  the  energy  norm  and  K  is  a  constant  determined  by 
Lh 

pQ  and  II  ^P II a, »  the  °°  norm  the  gradient  of  p(x,y).  Moreover,  the 
proof  clearly  indicates  why  one  should  expect  great  improvement  when  further 
"smoothing"  is  introduced. 

Thus  we  extend  the  results  of  Braes s  [2],  Ries,  Trottenberg  and  Winter  [8] 
to  include  a  variable  diffusion  coefficient  p(x,y)  and  more  general  regions. 

In  section  2  we  formulate  the  problem  and  the  basic  two-grid  method  of 
solution.  In  section  3  we  prove  the  basic  estimate.  This  proof  proceeds  from 
a  fundamental  insight  of  McCormick  and  Ruge  [7].  Section  4  describes  the  re¬ 
sults  of  some  computational  experiments  which  lead  one  to  believe  that  the 
results  of  [2.]  are  essentially  correct  for  the  variable  coefficient  case  as 
well.  These  computations  were  carried  out  on  the  CRAY  I  at  the  Los  Alamos 
National  Laboratory.  Finally,  an  appendix  gives  the  basic  "energy"  estimate 
required  in  section  3. 
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2.  The  Problem 

Given  a  (small)  value  h  >  0  let  {(x.,y.)  =  (kh,jh);  k,j  =  0, 

K  J 

+  1,  +  2,...}  be  the  associated  mesh  points  in  the  x  -  y  plane.  Let 

(2.1a)  RE:=  {(xk,yj);  k+j  =  0  (mod  2)}  , 

(2.1b)  RQ  :=  {(xk,yj);  (k+j)  =  1  (mod  2)}  . 

Let  SI  be  a  bounded  domain  in  the  plane  with  a  piecewise  smooth  boundary 
3ft  .  We  wish  to  define  the  set  of  "interior"  mesh  points,  ft^  .  We  assume 
that  h  is  less  than  j  the  length  of  each  smooth  section  of  3ft  . 

The  main  result,  (Theorem  3.1)  for  the  two-grid  iterative  schemes  is 
valid  in  quite  general  domains  provided  that  we  use  a  modification  of 
"approximation  of  degree  0"  (see  [6])  to  describe  the  boundary  conditions 

Definition:* 

(i)  If  (x|<>yj)  e  RQ  n  ft  we  say  that  (x^y^)  e  fth  if  the  four 

neighbors  ( (xk^.^  »yj ) »  (xk_i»yj)»  (^k»yj_,j)»  (xk,yj^_-j)} 

and  the  line  segments  from  (x. ,y.)  to  each  of  its  neighbors 

K  J 

all  lie  in  ft,  the  closure  of  ft. 

(ii)  If  (xk,y^)  e  R^  n  ft  we  say  that  (xk>yj)  e  ft^  if  the  eight 

neighbors  {(x^.y^.),  (xk_1,yj),  (x^y^),  (x^y^),  (xk+1,yj+1), 

(xk+ryj-l}’  (xk-l’yj+l*’  (xk-1  »yj_i )}  an<^  the  line  segments  from 
(xk,y. )  to  each  of  its  neighbors  all  lie  in  ft  . 


(*) 

'We  must  consider  the  line  segments  from  (xk,y.)  to  the  neighbors  only 
in  the  case  of  reentrant  corners  or  cusps.  K  J 


When  3 ft  has  a  cusp  or  a  corner  at  a  point  (x,y).  we  require  that 


(x,y)  =  (xk,y.)  e  re  . 

The  points  (x^.y^)  €  ft/ft^  are  the  boundary  points  of  ft^  .  That  is 

3fth  :=  { ( x k , y ^ )  £  Sl/Qh}  . 

A  true  multigrid  requires  the  use  of  many  coarser  grids.  In  such 
general  regions  the  treatment  of  the  boundary  conditions  on  succeeding 
coarser  grids  gets  complicated.  In  truth,  the  multigrid  literature  has 
barely  touched  on  this  question.  In  the  case  studied  by  Braess  [2],  ft 
is  a  polygonal  domain  whose  sides  have  slope  +1,0  or  °°  and  the 
corners  all  belong  to  the  coarsest  (and  hence,  the  finest)  grid.  For 
this  case  we  note  that  (see  Figure  1): 

(ftl.a)  3ft^  c  3ft 

and 

(ftl.b)  if  a  is  a  side  of  ft  with  slope  +  1,  then  all  the  points  of 
3ft^  which  also  lie  on  3ft  belong  to  . 

For  any  function  F(x,y)  defined  on  the  (x,y)  plane  we  write 


(2.2a) 

^k,j  :=  ^xk’^j^  ’ 

(2.2b) 

(2.2c) 

:=  F<V(jt,i>h)  ■ 

To  obtain  an  approximate  solution  of  (3)  we  seek  a  grid  function 
{Ukj}  defined  on  the  mesh  points  and  satisfying  the  system  of  equations 

f°r  (xk»yj)  e 


(2.3a) 


h2  {pk,j-*CUk,j'Uk,j-l] 


3k,j+J5CUk,j+TUkJ]} 


and,  for  (x^y^)  i 

(2.3b)  Ukj  =  0  . 

We  rewrite  (2.3)  as 

<2-4a>  Mu  ’  fkj  •  <V*J>  ‘  °h  • 


<2-4»)  Ukj  '  0  •  <V*J>  2  \  ' 

We  turn  to  the  question  of  the  solution  of  these  linear  algebraic 
equations  via  a  "two-grid"  method.  Let 

(2.5)  ft^  :=  R^  n  ft^  ,  ^lg  Rg  n  ^  • 

Our  two  grids  are  ft^  and  ft^  .  Let  Sh  and  be  the  spaces  of  grid 
functions  defined  on  Rr  u  R_  and  Rr  which  vanish  outside  ft.  and  ftr 
Our  first  step  is  to  set-up  "communication"  between  these  two  spaces. 

To  be  specific,  we  construct  linear  "interpolation"  and  "projection"  opera 
tors  1^  ,  so  that 

(2.6a)  1^:  Sh  -*•  ,  (Projection), 


(2.6b) 


’e:  SE  ~  V 


(Interpolation) . 
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Lemma  2.1 :  Let  lj|  be  defined  by  (2.7).  Then,  a  function  U  =  U(h)  e  Sh 
is  in  R  if  and  only  i f 

(2.10)  ^LhU\j  =  0  with  (xk,yj)  €  %  ‘ 

We  are  now  ready  to  describe  the  two-grid  methods.  Let  G  be  a  smoothing 
operator.  That  is,  given  u^  £  S.  we  construct  u  via 


(2.11a) 


u  =  Gu°  =  u°  +  B(f-Lhu°)  =  Gqu°  +  Bf 


(2.11b) 


G0  -  (I-BLh) 


where  B  is  a  given  matrix  and 


(2.11c) 


<LhGo^,G0^>  <  ! 


G_ ||  I"  *  sup  0^0y 
h  <*°  <Lh*,*> 


Algorithm  2.1  : 


Step  1 :  Given  u^  e  S^,  form  u  =  Gu^ 


Step  2:  Construct  u  via 


Skj*“kj'  6  re 

[Lh0]kj  =  fkj  •  (Vyj>  e  ao 

®kj  ’  0  ’  lxk’yj*  e  R0/t!0 


That  is:  "relax"  the  equations  on  the  "odd"  points, 


9 

* 

/v  E 

Step  3:  Set  r  =  f  -  L^u  ,  r^  =  Ihr  . 

Step  4:  Solve  L^<j>  =  r^  where  is  the  "coarse  grid  operator"  to  be 
I  described  later. 

1  ~  h 

Step  5:  Set  u  =  u  +  . 

i  Step  6:  Set  -*■  and  return  to  step  1. 

An  important  smoother  G  is  the  odd-even  Gauss-Seidel  scheme.  That 
is,  define  H°  -  relaxation  on  the  odd  points  as  in  Step  2  - 

i 

(2.12a)  {H°u)kj  =  ukj  »  {xk’yj^  e  RE  * 

and 

i 

(2.12b)  Cih(H°.nkJ  -  fkJ .  Oy^'V 

(2.12c)  uk>j.  =  0  ,  (xk,yj)  £  Rq/JIq  . 

Similarly  define  HE  ,  relaxation  on  the  even  points,  by 
(2.13a)  (HEu)kJ=bkj.  (Xk €  R0  ^  RE^^E  * 

i 

(2.13b)  CLh(HEu )]kj-  =  fkj  ,  (xk’yj)eJV 

Let  v  >  0  be  a  integer.  We  obtain  the  generalized  MGR[v]  two-grid  iterative 
scheme  by  choosing 

(2.14)  G  =  (HEH°)V  . 


We  now  describe  two  choices  of  the  coarse  grid  operator  LF. 


Case  1 :  Let 


(2  15a)  a  -  1  pk^,jpk-l,j-^  ,  pk,j-Japk-fr,J-1 

(  a)  h2  ck-l,j  ck,j-l 


(2.15b)  b.  ,  .  ,  =  X  Pk»J-^Pk+i5»j~1  +  Pk+fr.JPk+1  .J-?» 
k+**J‘*  7  ck,j-l  ck+l,j 


{2J5C)  dkj  =  [ak-Js>j-3s+ak-HsJ-U5+bk+J5,j-)5+bk-J5,j+l5;l 
Then,  if  (k+j)  =  0  (mod  2). 


(2.16) 


[LE  U]kj  =  '  ak+>s, j+*5Uk+l  ,j+l  "  ak-Js, j-^Uk-l  ,  j-1 


”  bk+>s,  j-%Uk+l ,  j-1  '  bk-%.j+JsUk-l ,  j+1  +  dkjUkj 


Case  2;  (The  Standard  Case):  if  k+j  h  0  (mod  2)  then 

( 2 )  1 

[LE  U]kj  =  Z1  {~pk+J5,  j+JsUk+l ,  j+1  "  pk+*s,  j-*sUk+l ,  j-1 
(2.17a)  2h 

~Pk-*s, j-*sUk-l , j-1  ’  Pk-Js, j+hUk-l  ,j+l  +  SkjUkj* 


where 


(2.17b)  Skj  =  {pk+is,j+J5+pk+is,j-J5+pk-J5,j+J5+pk-J5,j-J5}  ' 


\  -* •* .  *',"***  ’  •  -  *  •  i'-  **.  •*.  •*-  ■*.  **.  •* .  •*.  ’  »  *  •> 

•  ’  *  .*  *  •  *  * .■  '  •  S  V.  *  .  “  •  ‘  •  “  -  .  ‘  .  *  w  ♦  .  •  «  «*  .  ■  •  ‘  -  •  .  ■ 
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3.  Analysis  of  the  Algorithm 

We  begin  our  analysis  with  an  observation  which  is  (by  now)  well  known 
among  multigrid  theorists  (see  [7]).  Let 


(3.1) 


I^L  1^ 
h  h*E 


Consider  Steps  4-5  of  the  two-grid  iteration.  Suppose  we  replace  L£  by 
Le.  i.e.,  suppose  we  find  the  function  which  satisfies 


and  set 

1  ^  T  h , 

u  =  u  +  I^ij;  . 

We  claim  that 


i.e.  u1  is  the  desired  solution!  To  see  this  we  set 


(3.2a) 


/\ 

c 


U 


u 


and  observe  that  Step  2  implies  that  if  k  +  j  =  1 


(Lh£)kj  =  (LhU-Lhu^kj  =  (f'Lhu)kj 


Hence  Lemma  2.1  asserts  that  there  is  a  function 


(mod  2),  then 
=  0  . 

V  £  SE  and 


/s 

e 


(3.2b) 
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We  now  verify  that 


LeV  =  I^L^V)  =  !hLh^  =  rE  • 

Hence, 


<|>  -  V 


and 

(3.3)  u  -  I =  u  -  e  =  U ! ! 


Unfortunately  we  have  chosen  Step  4  with  and  not  .  This  choice  was 
not  merely  pique  on  our  part  (or  the  part  of  Braes  and  Ries,  Trottenberg  and 
Winter).  The  point. is  --  having  chosen  as  a  five  point  star  we  can  now 
proceed  to  replace  Step  4  with  a  new  two  grid  step  --  i.e.  we  can  build  a 
true  multigrid. 

In  any  case,  the  problem  of  Step  4  is  seen  to  be 

^  =  ’ 

where,  as  we  see  from  Lemma  2.1,  ij^  is  the  Lh  projection  of  e  into 
R  .  Hence 


(3.4a)  l|c|IL  -  ||  I^||L  <  ||S||,  . 

h  h  h 

We  will  give  a  complete  description  of  in  the  appendix.  For  now,  we 
wri  te 


where  is  defined  by  this  equation.  Observe  that  both  and 
(either  or  L^)  are  symmetric,  positive  definite  operators. 

Hence  the  associated  is  a  symmetric  operator.  Our  main  estimate  is 


Lemma  3.1:  For  '  or  LE  =  ,  there  is  a  constant 

only  on  II il <*,  *  the  maximum  norm  of  the  first  derivatives  of 
coefficient  p(x,y),  and  pQ  such  that,  for  all  <j>  e  SE>  <M 


<  Lp$ ,  4^ > 

(3.6)  -Kh<— - -  <  2(1+Kh)  . 

<  ) 

Proof:  See  Theorem  A  of  the  Appendix. 

Consider  the  eigenvalue  problem 
(3.7a)  (ALe  -Le)  <p  »  0  ,  tp  +  0 

which  is  equivalent  to 

(3.7b)  (AI  -  l£}lE)  ip  =  0  ,  ip  fO  . 

Using  (3.5)  we  see  that  this  problem  is  equivalent  to 
(3.8a)  [(2A-1)LE -LE]^  =  0  . 

From  Lemma  3.1  we  find 


(3.9) 


1-Kh  „  ,  ,  3+2Kh 

—  1  A  £  2 


Theorem  3.1:  Let 


c 


0 


then 


11^1,  <i(l+Kh)||c°||L  . 


K,  depending 
the  diffusion 
0  we  have 


(3.10) 


Proof:  We  have 


Using  (2.11c)  we  see  that 

(3.1D  ||  e||  L  il|e0||L- 

h  h 

From  (3.2a),  (3.3b)  and  Step  5  of  the  multigrid  algorithm  we  have 

(3.12a)  e1  =  e  -  I^4>  -  I^(ip-d>)  , 

and 

(3.12b)  Le4>  -  . 

Hence 

(3.13)  it  -  <f»  -  (I-L^Lg)^  . 

Thus 

(3.14a)  <  L^e1  .e1  >h  =  <  >h  =  2  <  I hLhrE  * ('J'-(f>)  \ 

=  2<  Le(^-<()),(^-4>)  >  E 

=  2<Le(I-L"1Le)^,(I.L^1Le)(/>e  , 

(3.14b)  ||  e1  ||  f  «  2<CI-L^Le1L^]L^.[I-^Le1L^]L^>e  . 

h 

_1  ~ 

Since  the  eigenvalues  of  LE  L  are  also  the  eigenvalues  of  the  symmetric 
operator  L^LE  ,  (3.9)  implies  that  the  eigenvalues  y  of  the  symmetric 

operator  ( I -L^L^ 1 L^)  satisfy 

-|(l+Kh)  <  y  <  |^(l+Kh)  . 


Thus,  (3.14b)  implies 


* 


1  n  2 


e'll[  1  2y2  <  L^,L^  >£  =  2y^<LE^,4»>E  =  2yi^<  Lhl”\p,(I^)  'ip>| 


rExT, 


=  2y2l<LhI^,IhE^>h 


=  U2<Lh£,£>h  =  y2||e||2 


This  result,  toget'er  with  (3.4a)  and  (3.11)  implies  the  Theorem. 


In  order  to  demonstrate  that  the  results  of  section  3  are  valid  for 


the  variable  coefficient  case  an  experimental  project  was  undertaken. 

The  essence  of  this  project  was  to  write  a  computer  program  which  imple¬ 
mented  Algorithm  2.1.  By  experimenting  with  different  functions  p(x,y) 
and  different  true  solutions  u(x,y)  it  was  shown  that  formula  (2)  of 
section  1  is  valid  for  the  variable  coefficient  case.  The  region  ft  is 
the  unit  square. 

The  computer  program  runs  in  an  interactive  fashion  and  allows  the 

user  to  provide  a  number  of  parameters.  These  include  N,  the  number 

of  points  on  a  side  of  ft^,  the  fine  grid  and  v,  the  number -of  smoothing 

iterations.  Starting  with  a  particular  initial  guess.  Algorithm  2.1  was 

then  repeated  until  the  discrete  norm  of  the  residual  was  less  than 
-8  0 

10”  .  For  the  initial  guess  U  ,  interior  points  of  ft^.  were  set  to  5 
while  interior  points  of  ft^  were  set  to  -5. 

Experiments  were  done  with  l^,  the  coarse  grid  operator  chosen  to  be 
both  and  .  The  calculation  of  was  complicated  by  the 

fact  that  for  points  of  ft^  for  which  refers  to  points  of  8fth, 

formula  (2.15c)  does  not  apply.  The  reason  for  this  is  because  the  computa¬ 


tion  of  either  or  b^±^._^  involves  referring  to  points  outside 

of  ft.  .  Of  course  since  U.  .  =  0  if  U.  .  e  8ft.  we  set  a...  ...  and 
h  kj  kj  h  k±*5j±J} 

to  zero  when  Uk±1j±1  and  Uk±lj.±1  are  in  8fth-  However,  we 
still  need  a  value  for  dkj.  for  the  two  nearest  interior  points.  For  the 
four  corners  points,  we  set  dkj.  to  be  the  value  of  dkj.  of  the  nearest 


interior  point.  As  the  mesh  gets  finer,  this  approximation  to  the  true 


d^  improves.  However,  in  almost  all  of  the  experiments  the  rate  of 

convergence  using  was  not  quite  as  good  as  the  rate  obtained 

(2) 

using  '  . 

The  tables  below  list  the  functions  p(x,y)  and  the  true  solutions 
u(x,y)  used  for  the  experiments.  For  each  problem  the  numerical  results 
obtained  using  both  and  are  displayed.  N  corresponds  to  the 

number  of  interior  points  on  a  side  of  and  v  corresponds  to  the  number 
of  smoothing  iterations.  The  smoother  used  was  the  odd-even  Gauss-Seidel 
scheme  as  described  in  section  2.  a(v)  in  the  tables  corresponds  to  the 
theoretical  rate  given  in  equation  (2)  of  section  1.  The  theoretical  rate 
has  c-.ly  been  proven  to  be  valid,  when  v  >  0,  in  the  constant  coef¬ 
ficient  case.  However,  as  can  be  seen  from  the  numerical  results  it 
appears  to  be  valid  in  the  variable  coefficient  case  as  well. 

In  conclusion,  the  numerical  results  demonstrate  the  validity  of 
Theorem  3.1  for  the  case  v  =  0  and  support  extending  equation  (2)  of 
section  1  to  the  variable  coefficient  case. 
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Table  I,  Experimental  Results 

Problem  1,  p(x,y)  =  1  ,  u(x,y)  =  0 

L  .L(,) 
lE  lE 


\  V 

0 

1 

2 

3 

15 

.4858 

.0646 

.0344 

.0200 

31 

.4844 

.0696 

.0375 

.0252 

63 

.4836 

.0708 

.0386 

.0263 

a(v) 

.5000 

.0741 

.0410 

.0283 

Problem  2,  p(x,y)  =  1  ,  u(x,y)  =  sin  ttx  sin  iry 
,  ..(1) 


X 

0 

1 

2 

3 

15 

.4858 

.0646 

.0344 

.0200 

31 

.4844 

.0696 

.0375 

.0252 

63 

.4836 

.0708 

.0386 

.0263 

a(v) 

.5000 

.0741 

.0410 

.0283 

L  =L{2) 

le  le 


BS 

0 

1 

2 

3 

15 

.4858 

.0646 

.0344 

.0200 

31 

.4844 

.0696 

.0375 

.0252 

63 

.4836 

.0708 

.0386 

.0263 

a(v) 

.5000 

.0741 

.0410 

.0283 

L  =L(2) 

le  e 


0 

1 

2 

3 

15 

.4858 

.0696 

.0344 

.0200 

31 

.4844 

.0696 

.0375 

.0252 

63 

.4836 

.0708 

.0386 

.0263 

a(v) 

.5000 

.0741 

.0410 

.0283 

a 


fr 


k; 

ia 


X  x  y 

Problem  6,  p(x,y)  =  e  (1  +  *s  sin  Try),  u(x,y)  =  e  J  sin  irx  sin  iry 


L  =  L(1) 
LE  le 


L  >L(2) 

le  le 


\  V 

0 

1 

2 

3 

15 

.4879 

.1084 

.0727 

.0565 

31 

.4854 

.0901 

.0582 

63 

.4851 

.0784 

.0473 

.0350 

0(v) 

.5000 

.0741 

.0410 

.0283 

\  V 

N\ 

0 

1 

2 

3 

15 

.4869 

.0686 

.0377 

.0255 

31 

.4851 

.0710 

.0390 

.0270 

63 

.4850 

.0715 

.0390 

.0270 

a(v) 

.5000 

.0741 

.0410 

.0283 

Problem  7,  p(x,y)  =  e~xy  ,  u(x,y)  =  (l-eA)(x-1)y  cos 


TTy 


L  =L(1) 
E  LE 


L  =L(2) 

le  le 


\  V 

0 

1 

2 

3 

15 

.4857 

.0797 

.0482 

.0351 

31 

.4842 

.0739 

.0431 

.0312 

63 

.4836 

.0714 

. 

.0399 

.0278 

a(v) 

.5000 

_ 

.0741 

_ 

.0410 

.0283 

Problem  8,  p(x,y)  =  e^sin  2  cos  ^  ,  u(x,y)  =  e"X'yx(x-l  )y(y-l ) 


BSI 

0 

— 

1 

2 

3 

15 

.4853 

.0650 

.0347 

.0207 

3' 

.4841 

.0697 

.0376 

.0253 

63 

.4835 

.0708 

.0386 

.0263 

a(v) 

.5000 

.0741 

-| 

.0410 

.0283 

Problem  9,  p(x,y)  = 


I  for  k  <  x  <  %  ,  74 -x  <  y  <  h  +  x  ,  x-J*  <  y  <  %-x 

6  otherwise 


u(x,y)  =  e“  yx(x-l )y(y-l ) 


\  V 

N 

0 

1 

2 

3 

15 

.4848 

.0648 

.0311 

.0174 

31 

.4857 

.0698 

.0377 

.0248 

63 

.4855 

.0711 

.0389 

.0266 

o(v) 

.5000 

.0741 

.0410 

.0283 

\  V 

0 

1 

2 

3 

15 

.4649 

.1815 

.1015 

.0634 

31 

. 

.4854 

.1186 

.1291 

.1046 

63 

- 

.4855 

.0712 

.0393 

.0275 

a(v) 

_ _ 

.5000 

.0741 

.0410 

.0283 

In  this  section  we  determine  and  the  quadratic  forms 

<  >  ,  <  L^ip.ip  >  ,  <  LEip,ip  > 

Let  U  e  Sc,  let  (x. ,y.)  e  Qr  .  Then 

t  K  J  t 

(A.D  Mkj-jtvKj- 

For  any  V  e  Sh>  [L^V]^.  involves  the  four  values  Vk±1  y  Vk  ^±1  . 
Therefore  we  consider  the  four  squares  I,  II,  III,  IV  (see  fig.  2)  with 
vertices 


(A. 2a) 

I: 

(xk+l’yj+l^’ 

(xk+2,yj)’  (xk+1 ,yj-l )}  • 

(A. 2b) 

II: 

l<Vyk1, 

(xk+l ,yj+l ^ 

(xk,yj+2)>  (xk-l  ,yj+l )}  ’ 

(A. 2c) 

III: 

{(xk’yj>' 

(xk-l ,yj+l 

(xk-2 ,yj ) ’  (xk-l ,yj-l ’ 

(A. 2d) 

IV: 

(xk-l ,yj-l 

(xk’yj-2)’  (xk+l ,yj-l • 

h 

In  each  square  the  value  of  [I^U]  at  the  center  point  is  a  weighted  average 

(given  by  (2.7b),  (2.7c))  of  the  values  of  U  at  the  corners.  Thus,  in  general 
Le  is  a  9-point  operator  based  on  the  9  vertices  of  these  four  squares. 

Since  LE  is  a  symmetric  operator  it  takes  the  form 


A 

[LEU]kj  =  EkjUkj  '  ak+l , jUk+2, j  ‘  ak-l,jUk-2,j 
■6k,j+lUk,j+2  "  Sk,j-lUk,j-2 

(A. 3) 

'Yk+Js,j+hUk+l  ,j+l  ‘  Yk-Js,  j-JsUk-l ,  j-1 


CTk-*5,  j+^k-l ,  j+1  '  ak+J5,j-J5Uk+l  ,j-l  . 


Lemma  A.l:  Let 


(A. 4a) 
where 


Ekj  =  ak+l ,  j  +  ak-l,j  +  gk,j-l  +  ^k, j+1  +  Yk+Ji,j+*5 

(A. 4b) 

+  Yk-*s,  j-J?  +  ak+Js,j-1i  +  ak-Js,j+i5  ‘ 

Then 


<Le^>  =  I  “k+1.j[^k+2,j -\jy 


(A. 5) 


+  ^  ^k,j+l  tyc,j+2  ”  ^kj^  +  ^  Yk+*5,j+ij^k+l ,  j+1 


+  ^  j-^k+l  ,j-l  ”*kj^  +  ^  Ekj^kj  * 


Proof:  Summation  by  parts. 


Similar  calculations  yield 


Lemma  A. 2:  Using  the  definitions  (2.15),  (2.17)  we  have 


<L^^,ij>>  =  E  ak-Us, j+^k+l , j+1  "  ^kj^ 


(A. 6) 
and 


+  ^  bk+i5, j.Jj^k+l  ,j-l  "^kj^  * 


<LE2)^’^>  =  ^7  Z  Pk+*s,  j+h^k+l , j+1  “*kj]2 


(A. 7) 


+  ?u2  ^  Pk+Js, j-Jj^k+l , j-1  '*kj]‘ 


We  now  compute  the  contribution  of  each  square  to  LL^lIJ ^ ^  and  the  two 
quadratic  forms.  In  evaluating  [L^I^U]^^  we  have  five  terms.  The  four 


terms 


'Pk±J5,j{IEU^k±l,j  *  'pk,j±Js{IEU)k,j±l 


are  clearly  associated  with  squares  (j[j)  and  (j*)  respectively. 


It  is 


convenient  to  agree  that 


pk±*s  jUkj  1s  associated  Wlth  square  (jjj)  . 
and 

pk  ji^kj  1s  associated  with  square  (jy)  . 

Let  Ekj(R),  ak±1 , j(R) *  8k,j±l(R)’  Yk±35,j±J5(R)’  °k±>s,jWR)  denote  the 
contributions  of  square  R  to  the  corresponding  coefficients  Ekj. ,  ak±1 

A 

6k,j±l’  Yk±Js,j±Js,  of  4  • 

Consider  square  I.  We  must  consider  two  cases,  either  (xk+1,yj)  e 
or  (xk+^,yj)  i  .  The  following  geometric  lemma  is  essential  to  under¬ 
standing  the  computations  in  the  latter  case. 


Lemma  A. 3:  Suppose  ( ,y ^ )  i  .  Then  either  ( xk+2 j )  t  n  or  the 
line  segment  from  (*k+1»yj)  to  (xk+2»yj)  is  not  entirely  in  fi  .  Further 

(xk+l,yj+l^’  ^xk+l’yj-l^  1  QE  ’ 

Proof:  (see  Figure  3).  Since  (x^y^)  e  the  points  (xk+1,yj+1),  (xk+1> 
(xk+l *yj-l )  €  ^  the  other  hand,  (xk+1 »y^ )  i  implies  that  either 

(xk+2»yj)  or  the  line  segment  from  Uk+1,y.j)  to  (xk+2»yj)  l  ft .  If 
(xk+2,yj^  i  ^  then  clearly  (xk+1,yj+1)  and  (xfc+1 ,yj_1 }  t  «E  . 
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If  (xk+2,yj)  e  Q  then  a  portion  of  3ft  crosses  the  line  segment 
xk+.j  <  x  <  x|<+2’  y  =  yj  •  ^  that  portion  of  the  boundary  continues  smoothly 

near  (xk+1,yj.),  then  the  line  segments  from  (x|<+1,yj.±1)  to  (xk+2»yj) 
are  not  entirely  in  ft.  Finally,  if  there  is  a  non-convex  corner  (x,y) 
near  (xk+]*yj)  that  corner  (x,y)  e  .  Hence  that  corner  must  be 
(xk+i »yj+i  )  or  (xk+i ,y.  which  is  therefore  not  in  ft^  .  The  other  one 
is  not  in  ft£  because  h  is  less  than  ^  the  length  of  smooth  segments 
of  3ft  . 

We  return  to  the  calculation  of  Ekj(I),  ak+1  j ( I ) ,  Yk+Jj  j+Jj(I),  and 

o i  ji  .  ,(I).  Square  I  does  not  contribute  to  the  other  coefficients. 

k+’i.J-v 

Case  1:  (xk+1.yj)  e  «h= 

A  straight  forward  calculation  yields 

(A. 8a) 


EyO)  ‘  -T 

J  2Ir 


k+*s,j 


Ck+1 ,  j 


(A. 8b) 


(A. 8c) 


k+1 ,  j 


Yk+*s,j+V 


k+Js.j-V 


1 

pk+5s,  j  pk+ 

3/z,j 

2h2 

Ck+1 , j 

1 

Pk+Jj, j  pk+l 

2h2 

Ck+1 ,j 

1 

Pk-H5,j  pk+l 

2h2 

Ck+1 ,  j 

(A.8d) 


(A. 10c) 


y  (m  -  i  nc+^j+i 

Tk*-W  ’  2h2  'k.j*, 


(A.IOd) 


k^-W  1  2h2  ck,j+l 


Case  2:  (xk.yj+1)  i  £y 


Using  arguments  similar  to  those  used  in  case  2  of  square  I  we  have 
(A. 11a) 


Ekj(H)  =  2h2  PkJ+Js  ' 


As  for  Yk+^  j+J^H)  and  CTk_*5  j 4.15 (II)  we  may  use  the  formulae  of  (A.  10c) 
and  (A.IOd).  Finally 


(A.Hb)  6kiJt,(1I)  =  EkJUI)  -  Y^^OI)  -  CTk-Jj. j+^t 1 1 )  • 


Because  l_E  is  symmetric  it  is  not  necessary  to  compute  the  contributions 
from  squares  III  and  IV.  We  now  make  a  similar  decomposition  of  the  coefficients 


of  .  Set 


(A. 12a) 

]  1 

ak+is,  j+i5(  1  5  =  ^2  ' 

(A. 12b) 

ak+*s,j+is(II)  =  ^2 

(A. 12c) 

bk+>5,j-15(I)  =  7  ' 

(A. 1 2d) 

bk-*5,j+i5(II)  =  ~2 

Ck+1 , j 


Pi-  ixl,  Pi 


k,  j+1 
Plrj.1.  i  Pux‘ 


k+1  ,j 
Pt,  Pi. 


't-  <4.1 


Let  <L^U,U>j  and  <L^  'U,U>j  denote  the  contribution  of  square  I 

A  (  1  ) 

to  the  quadratic  forms  <L^U,U>  and  <  L^  ;U,U>  respectively.  Then,  using 
lemma  A.l  we  see  that 


(A. 11a) 


f  <lem>i  ■  vi.j(I)hw,r*k.jl2  + 

Bk+l,j(I)^k+l,j+l  '  VlJ-l^  +  Yk+Js,j-H5{I)^k+l  ,j+l  ‘^k.j^ 
+  Yk+3/2,j-J5{I)^k+2,j '^k+lj-l^  +  ak+H,j-H(I)t4,k+l,j-l  ‘♦kj^ 
+  ak+3/2,j-Hs(I)^k+l,j+l  '  ^k+2, j^2  » 


and 


(A. lib)  U 


=  ak+Jj,  ^^k+1 ,  j+1  ’^kj^2 

ak+3/2 ’  j- 1  ^ C^k+2 ,  j  ‘^k+l.j-l^2  +  bk-Ms,  j-^s(I )C^k+l  ,j-l  -^k,j^2 
hk+VjW^k+l.j+l  ”^k+2,j^  • 


A  basic  inequality  is 
Lemma  A. 4:  Let  if/  «  S£  then 


(AJ2a)  Ok+1  f j+i 
(A. 12b)  [*k+1J*i 
(A. 12c)  C^k+1>j+1 


'^k+lj-l^2  -  2{^k+l,j+l  ~\+2,^2  +  ^k+2,j  ‘  Vl  ,j-1]  ] 
“*k+lfj-l^  -  2{[*k+lfj+l  ‘^k,^2  +  ^*k,j  “*k+l,j-l]  } 

’  ^k+1 , j-1  ^  -  ^k+l,j+l  '’W.jJ2  +  ^k+2,j '^k+l.j-l^2 


+  E'J'tJ.l  4  0.1 


Therefore  from  (A. 11a)  and  (A. 11b)  we  have 


<L<'W,  ■ak+l,j<n[*k+2,j-*k,j]‘ 


(A.  17) 


+  Bk+l,j(I)IA+l,j+l  "*k+l,J-1^ 


Thus,  we  have  established  the  left  hand  inequality  of  (A. 14).  Using  (A. 12), 


i.e.,  the  definitions  of  a|(+j5  bk+Jj  etc‘  and  (A-8b)»  the 


definition  of  ak+1  ^.(1)  and  (A. 10b),  the  definition  of  3k+^  j(I)  we  see 
that  there  is  a  constant  K  ,  depending  only  on  UVpjl^  and  pQ  such  that 


(A. 18a) 

(A. 18b) 
where 


j  <  i(l+Kh) 
aI  "  2 

<  i(l+Kh) 
aI  c 


(A. 18c)  aj  =  any  of  ak+3/2,j->5’  bk+V2,j+h}  ' 

Therefore  (A. 17)  and  (A. 12c),  (A. 13c)  and  (A. 17b)  yields  the  right  hand 
inequality  of  (A.14) . 

Corollary; 

(A. 19)  1  <Le<M>!  <  |(l+Kh)  . 

Lemma  A. 6:  Suppose  (x^  yj)  i  ^  .  Then  the  conclusion  of  lemma  A. 6  and 
its  corollary  hold. 


Proof:  Calculations  similar  to  those  of  lemma  A. 5  now  yield 


(A.20)  0  <  <  4*,*  >,  =  akt,  _ _ j  -  *k/  • 

Thus,  as  before,  the  left  hand  inequality  of  (A. 14)  holds.  In  this  case 
(A. 11a),  (A. lib)  and  (A. 10c),  (A.IOd)  imply 

<  (1+Kh)  . 
aI  “ 

Thus,  the  lemma  is  proven. 

Theorem  A:  There  is  a  constant  K,  depending  only  on  II V p || ^  and  pQ 
such  that,  for  both  choices  of  LE(=L^,L^)  and  associated  L£  ,  if 
<J>  t  0  and  <p  e  SE  we  have 

<  LP4>  ,4>  > 

_Kh  2(1+Kh)  • 

Proof:  The  arguments  which  give  lemma  A. 5  and  lemma  A. 6  extend  to  all  the 
squares,  II,  III  and  IV.  Thus,  those  lemmas  imply  that  the  theorem  holds 
for  =  L^.  The  case  of  follows  from  (3.5)  and  the  ob¬ 

servation  that 


(1-Kh)  (  ^ ill, ip  >  (1+Kh)  <  >  . 
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